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LAW OF LARGE NUMBERS FOR RANDOM DYNAMICAL 

SYSTEMS 

KATARZYNA HORBACZ AND MACIEJ SLECZKA 



Abstract. We cosider random dynamical systems with randomly chosen jumps. 
The choice of deterministic dynamical system and jumps depends on a posi- 
tion. We proove the existence of an exponentially attractive invariant measure 
and the strong law of large numbers. 



1. Introduction 
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'"^ ' In the present paper we are concerned with the problem of proving the law of 

rS^ ', large numbers (LLN) for random dynamical systems. 

j^ ' The question of establishing the LLN for an additive functional of a Markov 

process is one of the most fundamental in probability theory and there exists a rich 
literature on the subject, see e.g. the monograph of Meyn and Tweedie [17] and 
the citations therein. However, in most of the existing results, it is usually assumed 
that the process under consideration is stationary and its equilibrium state is stable 
in some sense, usually in the L^, or total variation norm. Our stability condition 
\^^ ' is formulated in a weaker metric than the total variation distance. 

OO ! The law of large numbers we study in this note was also considered in many 

^«0 ' papers. Our results are based on a version of the law of large numbers due to 

'!::j- I Shirikyan (see [19], [SO])- Recently Komorowski, Peszat and Szarek [12] obtained 

^D ' the weak law of large numbers for the passive tracer model in a compressible envi- 

^^ I ronment and Walczuk studied Markov processes with the transfer operator having 

spectral gap in the Wasserstein metric and proved the LLN in the non-stationary 
case [25]. 

K^ • Random dynamical systems [4], [6]take into consideration some very important 

5_j I and widely studied cases, namely dynamical systems generated by learning systems 

C^ ' [I] , [8] , [10] , [E] , Poisson driven stochastic differential equations [5] , [15] , [21] , [22] , 

iterated function systems with an infinite family of transformations [H] , [53] , [23] , 

random evolutions [3] , [18] and irreducible Markov systems [26] . 

A large class of applications of such models, both in physics and biology, is worth 
mentioning here: the shot noise, the photo conductive detectors, the growth of the 
size of structural populations, the motion of relativistic particles, both fermions 
and bosons (see [5], [TT], [13]), the generalized stochastic process introduced in the 
recent model of gene expression by Lipniacki et al. [J. 

A number of results have been obtained that claim an existence of an asymptot- 
ically stable, unique invariant measure for Markov processes generated by random 
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dynamical systems for which the state space need not be locahy compact. We con- 
sider random dynamical systems with randomly chosen jumps acting on a given 
Polish space {Y, g) . 

The aim of this paper is to study stochastic processes whose paths follow deter- 
ministic dynamics between random times, jump times, at which they change their 
position randomly. Hence, we analyse stochastic processes in which randomness 
appears at times to < ^i < ^2 < ■ • ■ We assume that a point xq ^ Y moves accord- 
ing to one of the dynamical systems Hi : M+ xY ^ Y from some set {Hi , . . . , IIjv}. 
The motion of the process is governed by the equation X{t) = Ili{t,xo) until the 
first jump time ii. Then we choose a transformation q^ : Y -^ Y from a family 
{qs : s G 5 = {1,...,A'}} and define xi — qs(ni{ti,xo)). The process restarts 
from that new point xi and continues as before. This gives the stochastic process 
{X{t)}t>o with jump times {^1,^2, • ■ •} and post jump positions {xi, X2, ■ ■ ■}■ The 
probability determining the frequency with which the dynamical systems 11^ are 
chosen is described by a matrix of probabilities [pij]^ -^j^, Pij '■ Y -^ [0,1]. The 
maps qs are randomly chosen with place dependent distribution. Given a Lipschitz 
function ip : X —i' R we define 

Our aim is to find conditions under which 5„('0) satisfies law of large numbers. 
Our results are based on an exponential convergence theorem due to Slgczka and 
Kapica (see (5) and a version of the law of large numbers due to Shirikyan (see 

[IS], HO]). 

2. Notation and basic definitions 

Let (X, d) be a Polish space, i.e. a complete and separable metric space and 
denote by Bx the cr-algebra of Borel subsets oi X. By Bb{X) we denote the space 
of bounded Borel- measurable functions equipped with the supremum norm, Cb{X) 
stands for the subspace of bounded continuous functions. Let Mfin{X) and Mi{X) 
be the sets of Borel measures on X such that /x(X) < 00 for /i e M.fin{X) and 
nix) = 1 for // £ M.i{X). The elements oi Mi{X) are called probability measures. 
The elements of AifiniX) for which ij,{X) < 1 are called subprobability measures. 
By suppfi we denote the support of the measure /x. We also define 

Mi{X) = {^Ie MiiX) : f L{x)fi{dx) < 00} 

Jx 

where i : X — !• [0, cx)) is an arbitrary Borel measurable function and 

Ml{X) = {Ai e Mi{X) : [ d{x,x)n{dx) < 00}, 

Jx 

where x G X is fixed. By the triangle inequality this family is independent of the 

choice of x. 

The space A^i(X) is equipped with the Fourtet-Mourier metric: 



ll/ii - mWpM = sup{| / f{x){p.i - ^l2){dx)\ : f e J"}, 
Jx 

lere 

T^{feCi,{X):\fix)-fiy)\<dix,y) and |/(x)| < 1 for x,y e X}. 
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Let P : Bb{X) — !> Bb{X) be a Markov operator, i.e. a linear operator satisfying 
Plx = Ix and Pf{x) > if / > 0. Denote by P* the the dual operator, i.e 
operator P* : Mfin{X) -^ Mfin{X) defined as follows 

P*fiiA) := [ PlA{x)fi{dx) for AeBx- 
Jx 

We say that a measure /i* G Ali(X) is invariant for P if 



Pf{x)fi^{dx) ~ / f{x)fi^{dx) for every f e Bb{X) 
X Jx 

or, alternatively, we have P* ^<t = /i,. An invariant measure /i is attractive if 



n— >oo 



lim / P"/(x) ^(dx) = / f{x)n{dx) for / G C,,(X), ^ £ Xi(X). 
"~" /x Jx 



By {Pj; : X e X} we denote a transition probability function for P, i.e. a family 
of measures Pj; G A^i(X) for x G X, such that the map x i-> 'Px{A) is measurable 
for every A G Bx and 

PI{x) - / ,/(y)Px(d?y) for X G X and / G Bb(X) 

Jx 

or equivalently P*/^(A) = Jy P2,(A)^(dx) for A G ;Sx and ^ G A^/i„(X). We say 
that a vector (pi, . . . ,pAr) where pi -.Y ^ [0, 1] is a probability vector if 

TV 

2^Pi{x) = 1 for .X G Y. 

i=l 

Analogously a matrix [pij\ij where pij : Y — )• [0,1] for i,j G {1,...,A^} is a 
probability matrix if 

N 

^Py(x) = l for xgF and i G {1, . . . ,iV}. 

Definition 2.1. A coupling for {Px : x G X} is a family {Bx,y : x, y G X} of 
probability measures on X x X such that for every B G Bx^ the map X^ 3 (x, y) ^-> 
'Bx,y{B) is measurable and 

Bx^y{AxX)=Px{A), Bx^y{XxA)=Py{A) 

for every x,y E X and A G Bx ■ 

In the following we assume that there exists a subcoupling for {P^; : x G X}, i.e. 
a family {Q^^y : x, y G A"} of subprobability measures on A^ such that the map 
(x,y) H^ Qj:^y(_B) is measurable for every Borel B C A^ and 

Q^,y(A X A) < Px{A) and Q:,,j,(A x A) < Py{A) 

for every x,y £ X and Borel A C A. 

Measures {Qa;.j, : x,y G A} allow us to construct a coupling for {P^, : x G A}. 
Define on A^ the family of measures {Rx,y : x, y G A} which on rectangles A x B 
are given by 

Rx,y{A xB) = ^-^—^{P^(A) - Cix,y{A X X)){Py{B) - Qx,y{X X B)), 
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when Qa;,y(X^) < 1 and Rx,j,(A x B) ^ otherwise. A simple computation shows 
that the family {Bx,y : x,y €i X} of measures on X^ defined by 

(2.1) B^^y ^ Clx,y + 'R.x,y for x,y £ X 

is a coupling for {P^; : x G X}. 

The following Theorem due to M. Sl§czka and R. Kapica (see [9]) will be used in 

the proof of Theorem l4.1l in Section 4. 



Theorem 2.1. Assume that a Markov operator P and transition probabilities 

{Qx,y : x,y £ X} satisfy 

AO P is a Feller operator, i.e. P{Cb{X)) C Cb{X). 

Al There exists a Lapunov function for P, i.e. continuous function L : X ^ 

[0, cxd) such that L is bounded on bounded sets, lim^^^oo L{x) = +oo and for some 

Xe (0,1), OO 

PL{x) < \L{x) + c for xeX. 
A2 There exist F C X^ and a G (0, 1) such that suppQx.y C F and 



(2.2) / d{u,v)Qx.yidu,dv) < ad{x,y) for {x,y) e F. 

A3 There exist S > 0, I > and v £ (0,1] such that 

l~QxAX^)<ld{x,yy 

and 

Qx,y{{{u,v) G X^ : d{u,v) < ad{x,y)}) > S 
for {x, y) e F 
A4 There exist /3 G (0, 1), (7 > and R> such that for 

k{ (xn, 2/„)nGNo ) == mf {^ G Nq : (a;„, y„) G F and i(a;„) + L{yn) < R} 

we have 

_ ~ 4c 

^x.yP " < C* whenever L{x) + L{y) < -, 

1 — A 

where Ej. y denotes here the expectation with respect to the chain starting from (x, y) 

and with trasition function {Ba;^^ : x,y € X}. 

Then operator P possesses a unique invariant measure /i* G A^f (X), which is 

attractive in A^i(X). Moreover, there exist g G (0, 1) and C > such that 

(2.3) \\P*^f, - M, Wfm < g"C(l + / L{x)fi{dx)) 

Jx 

for /I G M[{X) and n G N. 

We will also need a version of the strong law of large numbers due to A. Shirikyan 
([19). |20)). It is originally formulated for Markov chains on a Hilbert space, however 
analysis of the proof shows that it remains true for Polish spaces. 

Theorem 2.2. Let (51, J^, P) be a probability space and let X be a Polish space. 
Suppose that for a family of Markov chains {{X^)n>o,^x)x&x on X with Markov 
operator P : Bi,{X) — > Bh{X) there exists a unique invariant measure ^^, G A^i(A"), 
a continuous function v : X ^f R+ and a sequence (7n)nGN of positive numbers such 
that 7„ ^ as 7Z ^- cx) and 

\\P*"'5x - ^J.*\\FM < lnv{x) for x e X. 
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If 

oo 

C = y^ 7n < OO 



n=0 



and there exits a continuous Junction h : X ^f IR+ such that 

¥.^{v{Xl)) < h{x) for xeX,n>0, 

where Ej; is the expectation with respect to P^;, then for any x £ X and any hounded 
Lipschitz function f : X ^ R we have 



n-l 

lim -^/(X,-)= / f{y)^i.{dy) 



fc=0 

Pj; almost surely. 

3. Random Dynamical Systems 

Let (Y, g) be a Polish space, M+ = [0, +00) and / = {1, . . . , N}, S* = {1, . . . , K}, 
where N and K are given positive integers. 

Let Hi : M+ x y — >■ y, « £ /, be a finite sequence of semidynamical systems, i.e. 

Ili{0,x)=x for ieI,xeY 

and 

Ui{s + t,x) = Ui{s,{Ui{t,x)) for s,t eR+, i e I and X eY. 

We are given probabihty vectors pi -.Y ^ [0,1], i £ I,p^ : Y ^ [0,1], s G 5, a 
matrix of probabihties [pij]ij^i , pij : y — > [0, 1] , i,j € I and a family of continuous 
functions Qs : Y ^)- Y, s € S . In the sequel we denote the system by (11, q,p). 

Finally, let (il, S, P) be a probability space and {tn}n>o be an increasing sequence 
of random variables t„ : O ^^ R+ with to = and such that the increments 
Ai„ = tn — tn-i, n € N, are independent and have the same density g{t) = \e~^* , 
t> 0. 

The action of randomly chosen dynamical systems, with randomly chosen jumps, 
at random moments tk corresponding to the system (Il,q,p) can be roughly de- 
scribed as follows. 

We choose an initial point xo & Y and randomly select a transformation Hi from 
the set {Hi, . . . , 11^} in such a way that the probability of choosing Hi is equal to 
Pi{xo), and we define 

X{t) = ni{t,xo) for 0<t<ti. 

Next, at the random moment ti, at the point 0^(^1,2:0) we choose a jump qs from 
the set {qi, . . . ,qK} with probability Pg{Ili{ti, xq)). Then we define 

xi = qs(ni{ti,Xo)). 

After that we choose 11^^ with probability Pii-^{xi), define 

X{t) ^ni^{t-ti,xi) for ti<t<t2 

and at the point Hi^ (^2 — ti, xi) we choose qs^ with probability p^^ (Hi^ (t2 — ti, xi)). 
Then we define 

X2 = gsi(nij(t2 -ti,xi)). 
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Finally, given x„, n > 2, we choose 11^^ in such a way that the probability of 
choosing Ili^ is equal to Pi„_ii„(x„) and we define 

X{t) = Ili,^ {t - tn,Xn) for tn <t < tn+l- 

At the point Hi^{/S.tn+i,Xn) we choose q^^ with probability p^ (ni^(Ai„+i,a;„)). 
Then we define 

Xn+i = (7s„(ni„(Ai„+i,X„)). 

We obtain a piecewise-deterministic trajectory for {X{t)}t>a with jump times 
{ti,t2, ■ ■ .} and post jump locations {a;i,X2, . . .}. 

We may reformulate the above considerations as follows: Let {£,n}n>o ^^d {^n}n>i 
be sequences of random variables, ^„ : fi — ^ / and 7y„ : fi — >■ 5 and let {yn}n>i be 
auxiliary random variables, ?/„ : 17 — >■ y, such that 

¥{^0 ^ ilxo = x) =Pi{x), 

P(^„ = k\xn = X and ^„_i = i) == Ptfe(x), 

and 

"'ln = s\y„=y)=p^{y) 



for ri > 1, x, y e Y, k,i £ I and s e S* . 

Assume that {^n}n>o and {?/n}n>o are independent of {tn}n>o and that for every 
n G N the variables 771, ... , ?7„_i, ^1, . . . , ^„_i are also independent. 

Given an initial random variable ^0 the sequence of the random variables {xn}n>o, 

x„ : il — > F, is given by 

(3.3) Xn =qr,„{Tli„_^{tn-tn-l,Xn-l)) for 71=1,2,... 

and the stochastic process {X{t)}t>Q, X{t) : fi ^ F, is given by 

(3.4) X{t)=U^^_,{t-tn-l,Xn-l) for tn-l<t<t„, 71=1,2,... 

It is easy to see that {X{t)}t>o and {a;„}„>o are not Markov processes. In order 
to use the theory of Markov operators we must redefine the processes {X{t)}t>o 
and {Xn}n>o in such a way that the redefined processes become Markov. 

For this purpose, consider the space Y x I endowed with the metric d given by 

(3.5) d{{x,i),{y,j)) ^ Q{x,y) + gd{i,j) for x,yeY,i,jeI, 

where Qd is the discrete metric in /. 

Now define a stochastic process {^(i)}t>o, ^{t) : O — > /, by 

£,{t) = ^„_i for <„_i <t <tn, 7i = 1, 2, . . . 

Then the stochastic process {{X{t),^{t))}t>o, {X{t),^{t)) : fi ^ F x / has the 
required Markov property. 

In many applications we are mostly interested in values of the process X{t) at 
the switching points i„. Therefore, we will also study the stochastic discrete process 
(post jump locations) {(a;„,^„)}„>o , {x„,^n) : ^ ^Y x I. Clearly {(a;„,^„)}„>o 
is a Markov process too. 

We consider the stochastic process {(a;„,^„)}„>o , {xn,£,n) : il ^-Y x I, defined 
by p.iP " P.3p with the help of the system (11, 5,^1). We will need the following 
assumptions: 
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The transformations Hi : M+ x Y ^i- Y, i £ I and Qs : Y ^ Y, s £ S, are 
continuous and there exists x» £ 1" such that 

(3.6) / e^^^ g{qs{Iij {t, x^)), qsix^)) dt <oo for j£l, s £ S. 

The functions p^, s G S*, and pij, i,j G /, satisfy the fohowing conditions 
X! \Pij(^) - Ptjiy)\ < Lpg{x, y) for x,y <eY, i e I, 



(3.7) 



^\Ps{^) -Psiy)\< Lpg{x,y) for x,yeY, 



sGS 



where Lp, Lp > 0. 

We also assume that for the system {Il,q,p) there are three constants L > 1, 
a S M and Lq > such that 

(3.8) ^py(2/)p(n,(t,a;),n,(i,y)) <Le"*£)(x,2/) for x,y eY, i € I, t >Q 

3&I 

and 

(3.9) '^p^{x)g{qs{x),qs{y)) < Lqg{x,y) for a:,yGy. 
ses 

For a;, 2/ e y, i > we define 

/n(t,a;,j/) = {je/ : giU,it,x),U,it,y)) < Le'-'gix,y)} 

Iq{x,y) ^ {s e S : g{qs{x),qs{y)) < Lqg{x,y)} 

Assume that there are pq > 0,qo > such that : for every ii,i2 £ I, x,y £ Y 
and i > we have 

Y^ P'iiAx)P^2jiy)>Po, 

(3.11) 3^Iuit,.,y)^ 

Yl Ps{x)Ps{y) > qo- 

selq(x,y) 

Remark 3.1. The condition (|3.1ip is satisfied if there are io G /, so £ -S swc/i that 
g{U,^,{t,x),n,„{t,y))<Le°''g{x,y) for x,y€Y,t>Q, 
Qilso (x), qso (y)) < Lqg{x, y) for x,y &Y, 

and 

inf inf p^ig {x) > 0, 
(3.13) ^eixeY 

inf p„ (x) > 0. 

To begin our study of the stochastic process {(a;„,^„)}„>o consider the sequence 
of distributions 

7I„(A)=P((a;„,erOe^) for A e B{Y y. I), n>Q. 

It is easy to see that there exists a Markov-FeUer operator P : M. ^f M. such that 

7i„+i = PV-n for n > 0. 
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The operator P is given by the forraula 
(3.14) 

PKA)=Y.Y. / / \e~^'lA{qs{n,it,x)),j)p,,ix)p,{Ujit,x))dt^iidx,di) 

and its dual operator U by 

(•+00 
(3.15) Uf{x,i)^J2J2 / \e-^'f{qsiUjit,x)),j)p,,ix)p,{U,it,x))dt, 
jei sGS-^o 

where A is the intensity of the Poisson process which governs the increment At„ of 
the random variables {in}n>o- The operator P given by p.l4p is called a transition 
operator for this system. 

4. The main theorem 
Theorem 4.1. Assume that system {Il,p,q) satisfies conditions p.6p - p.lip . // 

(4.1) LLy + j<l. 

then 

(i) there exists a unique invariant measure /i, £ A4\{Y~kI) for the process (a;„, S,n)n>o, 

which is attractive in A4i{Y x I). 

(ii) there exist q G (0, 1) and C > such that for fi G A^}(y x /) and n G N 



\\P*"p - t^AlFM < g"C(l + / gix,x,)ti{dx)), 

where x, is given by liS.b]) . 

(Hi) the strong law of large numbers holds for the process {xn,Cn)n>a starting from 
(xq, ^0) (z Y X I , i.e. for every bounded Lipschitz function / : y x / — > M and every 
Xq (^Y and ^q (z I we have 

1 "~^ f 

lim -y2f{xk,^k) = f{x,Of^*{dx,d£,) 

"^xafy o.l'most surely. 

Proof of Theorem \4.1\ 
We are going to verify assumptions of Theorem 12.11 Set X — YxI,F — XxX 
and define 

Q(2;i,ii)(x2,i2)(^) = 
/•+00 

J2J2 / Xe~^'{p,,j{xi)p^{ll^{t,xi))Ap,,j{x2)Ps(nj{t,X2))}x 

xlA{{qs{lijit,Xi)),j),{q,{nj{t,X2)),j))dt 
for A C X X X , where a Ab stands for the minimum of a and b. 

AO. The continuity of functions pij,Pg,qs implies that the operator P defined 
in p.l4p is a Feller operator. 



LAW OF LARGE NUMBERS FOR RANDOM DYNAMICAL SYSTEMS 

Al. Define L{x,i) — gi(a;,a;*) for {x,i) G X. By p.l5|) we have 
UL{x,i) <J2Y1 / Q{qs{n,it,x)),q.s{U,it,x,)))\e-^'p,,ix)p,{U,{t,x)) dt 

+ X!X! / Qiqs{nj{t,x^)),qsix^))Xe~^*Pt.iix)p^{Uj{t,x)) dt 

+ X!X! / Qiqsix*),x^)Xe~^*p^j{x)p^(Uj{t,x))dt. 
jei ses-^o 



Further, usmg p.6p . p.Sp and p.9p we obtain 

(4.2) UL{x,i) <aL{x,i) + b, 

where 

XLL„ 
a — 



A — a' 

(4.3) ^+00 

& = ^^ / Ae"^*g(gs(nj(i,x,)),gs(x*))dt + ^g(gs(x*),a;,), 

so L is a Lapunov function for P. 

A2. Observe that by ((X5|) . dSH]) and ([53]) we have for (a;i,ii), (^2,^2) e ^ 

C?(w,f)Q(a;i,ii)(2:2,i2)(du,dw) = 

Zm / Xe''^'{p,,^{xi)p,(Uj{t,xi)) Ap,,j{x2)Ps(lljit,X2))}x 

X e(g5(nj(t,xi)),g5(nj(t,a;2)))dt 

^+00 
^J2Y. \e-^'p^,Jixl)p,{U,it,Xl))Qiq,{U,it,Xl)),qs{U,{t,X2)))dt 

< l3g{xi,X2) < I3d{{xi,ii), (x2,i2)) 
with /3 = ^^ < 1 by Kn . 

A3. From (Ejl) and ^M it follows that 



1 -^^{pi,j(a;i)p^(nj(i,xi)) Api^j{x2)Ps{'nj{t,X2))} 
jei s<£S 

jei ses 

< ^^Kij(a;i)|p^(nj(t,a;i)) ~p^{Uj{t,X2))\ 
jei ses 

+ ^^Ps(nj(i,a;2))biu(2:i)-p»2j(a;2)| 

< LLpe"*g{xi,X2) + Lpg{xi,X2) +2Ngdiii,i2) 
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and consequently 

1-Q(xi,»i)(^2.»2)(^')<(ip + 1 -)g{xi,X2) + 2Ngd{ii,i2)- 

A — (y 

Fix xi,X2 £Y and ii,i2 G /• Define B = {((u, j), {v,j)) : g{u,v) < l3g{xi,X2),j S 
/}. If a > then there exists Tq > such that LL,e°* < ^ iov t < Tq. Set 
A = (0,To). If a < then there exists Tq > such that LL,e"* < (3 ior t > Tq. 
Set A = (TojCX)). In both cases define r = J^Xe^^*^ dt. For all a;,y G F, i G A, 
j e In{t,x,y) and s G /g(nj(i,2;),nj(t, j;)) we have 

(4.4) {{qs{n,{t,x)),j),{qs{n,{t,y)),j)) e B. 

From p. lip and (|i^ we obtain 

> I \e-^' E E feu(a^i)P,(nj(i,xi))Ap,,j(x2)p,(n,(t,a;2))}x 
xlB((g,(nj(t,a;i)),j),(g,(n,(t,a;2),i)))di 

^ I \e-^' Y, J2 {p^,J{Xl)p,{U,it,Xl))Ap,,,{x2)Ps{Il,{t,X2))}dt 

- / -^^"^* XI 1^ feij(a;i)K2j(a;2)Ps(nj(i,xi))p^(nj(i,a;2))}di 

> Poqor > 0, 



^ ]eIn(t,xuX2) se7,(nj(t,a:),nj(t,t/)) 



so A3 is satisfied. Since F = X x X, assumption A4 is trivially satisfied. 
From Theorem [2?l1 we obtain (i) and (ii). Set v{x,i) — C{g{x,x^,) + 1) and h{x,i) — 
C{g{x,x^) + 1 + j4^) for a; G X, i e /, with a,b as in (|4.3p . Iterating (|4.2p we 
obtain 

^XQ,oi'"i^n,£,n)) <h{xo,^Q) for a;o ^ X, ^q £ ^. 
Application of Theorem 12.21 ends the proof. 



The next result describing the asymptotic behavior of the process (a;n)n>o on Y 
is an obvious consequence of Theorem l4.1l Let jlo be the distribution of the initial 
random vector xq and /i„ the distribution of x„, i.e. 



fln{A) = P(x„ G A) for A G BY,n > 1. 



Theorem 4.2. Under the hypotheses of Theorem \4.1\ the following statements hold: 
(i) there exists a measure /2* G M\{Y) such that for any fiQ the sequence {fin)n>o 
converges weakly to /2,. Moreover, if 

P(xo eA)= fi^A) for AeBy 

then jlniA) = p.*{A) for A G By and n > I. 
(ii) there exist q G (0, 1) and C > such that 

Wfj'n- P-*\\fm <q'''C{l+ / g{x,x^)jla(dx)) 



for any initial distribution JIq G M\{Y) and n>l. 

(Hi) for any starting point xq ^ Y , S,o & I and any hounded Lipschitz function f on 

Y 



1 f 

lim -y^f{xk)^ / f{x)jl^{d2. 
fc=0 •'^ 



r(j,jQ almost surely. 
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[12: 
[is; 

[14 

[is; 
[16; 
[17: 

[is; 
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